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Here we report the combined spacecraft observations of Saturn acquired over one Saturnian year (—29.5 
Earth years), from the Voyager encounters (1980-81) to the new Cassini reconnaissance (2009-10). The 
combined observations reveal a strong temporal increase of tropic temperature (—10 Kelvins) around the 
tropopause of Saturn (i.e., 50 mbar), which is stronger than the seasonal variability (— a few Kelvins). We 
also provide the first estimate of the zonal winds at 750 mbar, which is close to the zonal winds at 2000 mbar. 
The quasi- consistency of zonal winds between these two levels provides observational support to a 
numerical suggestion inferring that the zonal winds at pressures greater than 500 mbar do not vary 
significantly with depth. Furthermore, the temporal variation of zonal winds decreases its magnitude with 
depth, implying that the relatively deep zonal winds are stable with time. 

T he Pioneer and Voyager spacecrafts made flyby visits to the giant planets and conducted snapshot observa- 
tions. The Galileo and Cassini spacecrafts on the other hand made orbital observations of Jupiter and Saturn 
respectively, which can be used to track seasonal variations. However, no single spacecraft observation is 
long enough to explore the temporal variation beyond the seasonal scale on the giant planets due to their long 
orbital periods. The Cassini observations in 2009-10 and the Voyager observations in 1980-81 both occurred in 
the northern early spring but were separated by one Saturnian year. Therefore, the combination of observations 
from Voyager and Cassini makes it possible to examine the non- seasonal variability of the atmosphere over one 
complete Saturnian year (—29.5 Earth years). Such an exploration will provide important information on the 
temporal variation of Saturn’s atmosphere, which is unique and different from the solar-forced seasonal variation. 

Here, we focus on exploring Saturn’s atmosphere in the tropical region (i.e., 30°N-30°S), in which the strong 
atmospheric winds —400-500 m/s blow from west to east. Our exploration of atmospheric temperature is mainly 
based on the thermal observations from the infrared spectrometer on Voyager and Cassini 1,2 . In addition, some 
independent radio -occultation measurements 3-6 are used to valid the observations from the infrared spectro- 
meter. Finally, the zonal winds measured by the imaging system on Voyager and Cassini 7,8 are used as the 
boundary condition in a modified thermal wind equation 9-11 to explore the vertical structure of zonal winds in 
Saturn’s upper troposphere (i.e., 50-750 mbar). 

In order to explore the temporal variation of Saturn’s atmosphere from one Saturnian year to the next, we 
searched all Cassini observations to select the best data and those as close as possible to one Saturnian year after 
the Voyager encounter to rule out the seasonal variation on Saturn. The selected observations are listed in Tables 
SI -S3 in the Supplementary Information. The data processing and the corresponding measurements for these 
selected observations are also discussed in the Supplementary Information. 

Results 

Figure 1 displays the tropic temperature in the upper troposphere from the nadir observations recorded by 
Voyager and Cassini. The zonal-mean temperature from the Infrared Spectrometer (IRIS) on the Voyager 
(1980-81) and the Composite Infrared Spectrometer (CIRS) on Cassini (2009-10) are displayed in panel A 
and panel B, respectively. Both of the Voyager/IRIS and Cassini/CIRS temperatures are retrieved from the 
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Figure 1 | Temporal variation of Saturn’s temperature from the Voyager 
epoch to the Cassini epoch. (A) Temperature retrieved from the Voyager/ 
IRIS nadir observations (1980-81). (B) Temperature retrieved from the 
Cassini/CIRS nadir observations (2009-10). (C) Differences of Saturn’s 
temperature between Voyager and Cassini. 


spectral measurements within the S(0) and S(l) hydrogen absorption 
lines with the same retrieval algorithm developed in a previous 
study 12 . The retrieval errors, which are mainly determined by the 
number of spectra, the signal strength, and the observational 


geometry, are estimated as — 1 Kelvin for the retrieved temperatures 
of both Voyager and Cassini 12,13 . The other errors due to the uncer- 
tain optical depth and He/H 2 ratio will probably introduce —2- 
Kelvin errors in the retrieved temperature (Section SI in the 
Supplementary Information). 

Panel C of Fig. 1 is the temporal variation of atmospheric temper- 
ature from the Voyager epoch to the Cassini epoch. This panel shows 
that there is significant tropical warming with a magnitude of — 10 
Kelvins around the tropopause from one Saturnian year to the next, 
which is significantly larger than the errors. Below the tropopause, 
the atmosphere shows a relatively small temporal variation. There is 
warming of —2-4 Kelvins around 15° in each hemisphere at the 
pressure level —300 mbar, which suggests that the temperature 
knee 14 ’ 15 probably changed its magnitude from the Voyager epoch 
to the Cassini epoch. At the pressure levels of 150 mbar and 
750 mbar, there is cooling of —2 Kelvins around the equator (i.e., 
10°S-10°N), which is close to the retrieval errors. 

Figure 2 displays the independent radio- occultation measure- 
ments in the epochs of Voyager and Cassini 5,6 , which are used to 
validate the temporal variation of atmospheric temperature revealed 
by the nadir observations from the infrared spectrometers on 
Voyager and Cassini (Fig. 1). The procedure of computing the tem- 
perature profiles from the radio-occultation measurements and the 
corresponding error estimates are introduced in Section S2 in the 
Supplementary Information. Figure 2 shows basic consistency of the 
atmospheric temperature around the tropopause (i.e., 50-70 mbar) 
between the nadir observations and the occultation measurements, 
which provides more evidences for the temporal variation of atmo- 
spheric temperature revealed by the nadir observations (Fig. 1). 
There are some discrepancies between the nadir observations and 
the occultation measurements existing below the pressure level of 
70 mbar, but most the discrepancies are smaller than the errors of the 
corresponding measurements. However, in the Voyager epoch the 
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Figure 2 | Comparison between the radio-occultation measurements and the nadir observations in the epochs of Voyager and Cassini. 

(A) Temperature comparison at 3.0°S amongthe radio-occultation measurements by Voyager 1 (1980), the nadir observations by Voyager 1/2 (1980-81), 
and the nadir observations by Cassini (2009-10). (B) Temperature comparison at 31.2°S among the radio-occultation measurements by Voyager 2 
(1981), the nadir observations by Voyager 1/2 (1980-81), and the nadir observations by Cassini (2009-10). (C) Temperature comparison at 1.8°S among 
the radio-occultation measurements by Cassini (2010), the nadir observations by Voyager 1/2 (1980-81), and the nadir observations by Cassini 
(2009-10). (D) Temperature comparison at 2.6°N among the radio-occultation measurements by Cassini (2010), the nadir observations by Voyager 1/2 
(1980-81), and the nadir observations by Cassini (2009-10). The latitudes shown in this figure are plenatographic latitudes. 
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discrepancies are larger than the corresponding errors in the pressure 
levels below 100 mbar for the latitudes of 3.0°S (panel A) and 31.2°S 
(panel B), which are probably due to some unknown systematic 
errors or underestimated errors. As we discussed in the Supplemen- 
tary Information (SI), most of the errors in the retrieved temperature 
are systematic errors. These systematic errors will affect the absolute 
magnitude of the retrieved temperature, but they do not affect the 
temperature gradients used in the thermal wind equation and hence 
do not alter our computations of thermal winds (please refer to S4 in 
the Supplementary Information). 

The significant temporal variation of atmospheric temperature 
over one Saturnian year is a new discovery. The maximal temporal 
variation of 10 Kelvins, which is larger than the seasonal variation of 
a few Kelvins in the same region 1316 , suggests that there are strong 
temporal changes besides the seasonal variation on Saturn. The tem- 
poral variation of tropospheric temperature is also different from the 
semi-annual oscillations 17,18 , because the semi-annual oscillations 
were discovered in the middle and high stratosphere of Saturn. 

The temperature fields retrieved from the nadir observations 
(Fig. 1) are further applied to a thermal wind equation to explore 
the vertical structure of zonal winds, as discussed in some previous 
studies 16,19,20 . Here, we used a modified thermal wind equation 911 
instead of the classical thermal wind equation 18-20 to explore the 
vertical structure of zonal winds in the tropic region of Saturn with 
emphasizing the temporal variation from the Voyager epoch to the 
Cassini epoch. In order to integrate the temperature with the modi- 
fied thermal equation, a boundary level with the known zonal winds 
is needed. Here, the pressure level of visible clouds at —360 mbar, in 
which the zonal winds are measured by the Imaging Science 
Subsystem (ISS) on Voyager and Cassini 7,8 (Fig. SI in the Supple- 
mentary Information), is used as the boundary level. 

Figure 3 shows a comparison of zonal -wind vertical structure 
between the Voyager epoch and the Cassini epoch. The errors in 
the computed zonal winds, which are —5-30 m/s, are estimated 
from the uncertainties in the cloud-tracking zonal winds at the 
boundary level, the uncertainties in the altitude location of the 
cloud-tracking zonal winds, and the propagating temperature errors 
in the integration of the modified thermal wind equation (Section S4 
in the Supplementary Information). Panel C of Fig. 3 suggests that 
the zonal winds around 10° latitude in each hemisphere decreased 
— 100 m/s at 50 mbar since the Voyager epoch. The — 100 m/s wind 
variation tilts equatorward from 10° latitude at 50 mbar to 5° latitude 
at 750 mbar, which suggests that the strong decrease of zonal winds 
discovered at the cloud level 21 basically exists in the whole upper 
troposphere (i.e., 50-750 mbar). 

Panel C of Fig. 3 also suggests that the temporal variation of tropic 
winds decreases its magnitude from 50 mbar to 750 mbar. The tem- 
poral variation of zonal winds at different pressure levels is further 
explored in Fig. 4. Panel A is the comparison of zonal winds at 
750 mbar, which shows that the magnitude of the equatorial jet 
(5°S-5°N) in the Cassini epoch (—450 m/s) is close to that of the 
equatorial jet in the Voyager epoch (—460 m/s). Such a comparison 
suggests that the 750-mbar zonal winds are relatively stable with 
time. Panel B of Fig. 4 further displays a comparison of the temporal 
variation of zonal winds between 360 mbar and 750 mbar, which 
shows that the temporal variation of zonal winds is stronger at 
360 mbar than at 750 mbar. Therefore, the relatively time-constant 
zonal winds at 750 mbar imply that the atmosphere with pressures 
greater than 750 mbar probably has relatively small temporal 
variability. 

Besides the temporal stability, the zonal winds at the relatively 
deep pressure levels also display stability in the vertical direction. 
The speed of the equatorial jet at 750 mbar (—450-460 m/s) is close 
to that of the equatorial jet at 2000 mbar (—500 m/s) 22,23 , which 
suggests that the equatorial jet does not vary significantly at pressures 
greater than 750 mbar. The nearly constant equatorial jet from 
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Figure 3 | Temporal variation of zonal winds from the Voyager epoch to 
the Cassini epoch. The zonal winds in Saturn’s upper troposphere (i.e., 
50-750 mbar) are computed from the temperature fields by a modified 
thermal wind equation. The regions within the black lines are left blank 
because the cylindrical routines of the modified thermal wind equation do 
not pass the regions. (A) Zonal winds derived from the Voyager/IRIS 
temperature (1980-81). (B) Zonal winds derived from the Cassini/CIRS 
temperature (2009-10). (C) Differences of Saturn’s zonal winds between 
Voyager and Cassini. 

750 mbar to 2000 mbar provides an observational evidence for a 
numerical inference on Saturn 24 in which zonal winds are suggested 
to extend from 0.5 bar to — 10 bar without significant decay, similar 
to their observed behavior on Jupiter 25 . 

Discussion 

In this study, the combined observations from Voyager and Cassini 
reveal that the tropical atmospheric temperature around the tro- 
popause warmed more than 10 Kelvins from one Saturnian year 
(1980-81) to the next (2009-10), which paints a remarkable picture 
of temporal variation of the thermal fields on Saturn. As the bound- 
ary of troposphere and stratosphere, the tropopause plays an import- 
ant role in planetary atmospheres. The significant warming in the 
tropopause will modify the atmospheric stratification and stability, 
and further influence the large-scale dynamics in the upper tro- 
poshere. Based on the temperature fields, we also explore the related 
wind fields, which provide some important information about the 
temporal stability and vertical structure of zonal winds in the upper 
troposphere. These important characteristics of the large-scale tem- 
perature and wind fields need to be considered in the future studies of 
the atmospheric system on Saturn. 

The 1990 tropical storm 26 is one possible driver for the large-scale 
atmospheric variation on Saturn. But the follow-up studies suggest 
that such a storm is probably not related to the significant warming 
around the tropopause (Fig. 1). First, the 1990 giant storm erupted 
around 5°N 26 , which is not consistent with the latitudes (20 °S and 
30°N) where the maximal warming happened. Second, a model 
study 27 suggests that the 1990 storm mainly developed below the 
tropopause, which probably did not affect the strong warming 
around the tropopause. Finally, numerical simulations 28 imply that 
the 1990 storm could not result in the large wind decrease — 100 m/s, 
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Figure 4 | Comparison of zonal winds between different times and between different pressure levels. (A) Comparison of750-mbar zonal winds between 
the Voyager epoch (1980-81) and the Cassini epoch (2009-10). The 750-mbar zonal winds come from the thermal winds shown in Fig. 3. (B) Comparison 
of temporal variation in the zonal winds between 3 60 mbarand750 mbar. The Temporal variation of zonal winds at 360 mbarand750 mbar comes from 
Fig. SI and Fig. 3, respectively. There are no Voyager observations for the latitude range of 2.0°S-10°S, which is interpolated from the neighboring 
measurements. Therefore, there are no error estimates for such a latitude range (2.0°S-10°S) (Section S4 in the Supplementary Information). The gaps 
around the equator (3°N-3°S) in the 750-mbar profiles are due to the integration routines of the modified thermal wind equation. 


which is related to the temperature changes. Even though the 1990 
storm is not directly related to the temporal variation discussed in 
this paper, it is possible that some dynamical changes induced by the 
1990 storm and other storms play roles in the large-scale atmo- 
spheric variation from Voyager to Cassini. 

In addition, the wave activity is the other possible driver for the 
large-scale atmospheric variation on Saturn. Waves are active on 
Saturn 29 . In the tropical region, we just identified a planetary- scale 
wave in the stratosphere 11 . In addition, the thermal maps from the 
Cassini/CIRS show some wavelike structures in the tropic region of 
troposphere (Fig. S4 in the Supplementary Information). Atmos- 
pheric waves can generate the inter- annual variability in the wind 
and temperature fields in the stratospheres of Earth 30 and Jupiter 31 . It 
is possible that the waves also play an important role in driving the 
large-scale variation in Saturn’s atmosphere from one year to the 
next. 

Finally, the energy budget affects the thermal structure on planets. 
Therefore, the exploration of energy budget provides one more per- 
spective to explore the temporal variation of the thermal structure 
and the corresponding dynamical fields. Saturn’s energy budget is 
mainly determined by the absorbed solar energy, emitted thermal 
energy, and internal heat. The time scale of the internal heat on the 
giant planets is generally thought to be much longer than a Saturnian 
year 32 , but the other two energy components (i.e., the absorbed solar 
energy and the emitted thermal energy) are variable on the time scale 
of Saturnian year. Our related study 33 revealed that Saturn’s emitted 
power varied from the Voyager epoch to the Cassini epoch. In addi- 
tion, Saturn’s absorbed solar energy probably changes with time 
because of the varying total solar flux 34 and the varying albedo due 
to the temporal variation of clouds 35,36 . We are exploring the tem- 
poral variation of the absorbed solar energy and hence the temporal 
variation of the energy budget on Saturn. Such a study will help us 
examine the role of energy budget in the temporal variation of atmo- 
spheric temperature revealed in this study. 


Methods 

Our results are based on the measurements of zonal winds and atmospheric tem- 
perature on Saturn in the two epochs (i.e., Voyager (1980-81) and Cassini (2009- 
10)). In addition, a modified thermal wind equation 9 ' 11 is applied to explore the 
vertical structure of zonal winds. In this section, we introduce the methods of mea- 
suring and computing the zonal winds and atmospheric temperature on Saturn. The 
full details of the data sets and the corresponding methods are included in the 
Supplementary Information. 

The retrieval of Saturn’s atmospheric temperature (Fig. 1) is mainly based on the 
infrared spectra recorded by the Infrared Spectrometer (IRIS) on Voyager 1 and the 
Composite Infrared Spectrometer (CIRS) on Cassini 2 . The process of retrieving 
atmospheric temperature is based on a retrieval algorithm developed by our co- 
authors and described in our previous study 12 . We also put some introduction of the 
retrieval algorithm in the Supplementary Information (Section SI). The retrieval 
errors of the retrieved temperature are estimated by considering these dominant 
factors affecting the retrieval, which include signal strength, varying optical depth, 
and the helium volume mixing ratio (please refer to Section SI in the Supplementary 
Information for the detailed discussion). To validate the retrieved temperature from 
the infrared spectrometers, we also use some independent observations from the 
radio-occultation measurements (Fig. 2) by Voyager (1980-81) and Cassini (2010), 
which come from previous analyses 5,6 . The detailed information of the previous 
analysis and the corresponding error estimates are put in Section S2 in the 
Supplementary Information. 

The measurements of Saturn’s zonal winds (Fig. SI) are based on the observations 
from the Imaging Science Subsystem (ISS) on Voyager 7 and Cassini 8 , which are 
described in Section 3 of the Supplementary Information. The basic method of 
measuring zonal winds is to track the visible cloud features on the ISS images. The 
errors of wind measurements are estimated by the standard deviation of the wind 
measurements in each latitude bin. Our discussion in Section 3 of the Supplementary 
Information shows that the other possible sources of errors in the wind measure- 
ments, which include the errors due to the navigation procedure and the uncertainty 
related to locating cloud features when conducting the cloud-tracking wind mea- 
surements, are negligible to the standard deviation of the wind measurements. 

The vertical structure of Saturn’s zonal winds (Fig. 3) are based on the integration 
of our modified thermal wind equation, which was described and validated in our 
previous studies 9 ' 11 . The basic methodology is to integrate the thermal wind equation 
with the retrieved temperature from the Voyager/IRIS and the Cassini/CIRS (Fig. 1) 
from a pressure level with known zonal winds. Such a pressure level is set as the 
pressure level in which the cloud-track zonal winds are measured from the ISS images 
(Section 3 of the Supplementary Information). The errors in the computed zonal 
winds (Fig. S5) from our modified thermal wind equation are estimated by the 
combination of three dominant error sources, which include the uncertainties in the 
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cloud-tracking zonal winds (Fig. S2), the uncertainties in the altitude location of the 
cloud-tracking winds (Figs. S3), and the uncertainty related to the errors in the 
retrieved temperature (Fig. S4). We discussed the three dominant error sources in 
Section 4 of the Supplementary Information, and the corresponding error estimates 
are used in Fig. 4 in the text. 
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